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The interpretation of Drell-Yan production in terms of the antiquark densities 
depends on NLO corrections. Besides the NLO corrections to the familiar annihi- 
lation — > 7* — » , there is a substantial contribution from the QCD Compton 
subprocesses gq — > 57* — > ql^l~ and gq — > 57* — > ql^l~ . The beam and target 
dependence of the two classes of corrections is different. We discuss the impact of 
this difference on the determination of the d — u asymmetry in the proton from 
the comparison of the pp and pn Drell-Yan production. 



The substantial d — u asymmetry of the nucleon sea is a striking manifesta- 
tion of the leading twist nonperturbative-tnesonic effects in parton distributions 
(for the recent review and references seel^l) . The first indirect experimental ev- 
idence |£or this asymmetry has been deduced from the Gottfried sum rule eval- 
uationH, but it left of course open the x— dependence of the d — u asymmetry. 
Following a suggestion by Ellis and Stirling , the CERN NA51 collaboration 
has measured the d — u asymmetry at a; c^i 0.18 from the comparison of the 
pp and pn{pd) Drell-Yan production. The asymmetry proved to be large and 
plaxjed at large x, in good agreement with expectations from meson cloud mod- 
els BQ. Recently, a much more detailed and higher precision comparison of the 
pp and pd Drell-Yan production has been performed in the Fermilab E866 
experiment. Furthermore there are plans of the dedicated Drell-Yan experi- 
ment at Fermilab with the goal of measuring the flavour content of the sea up 
to a; 0.(H. The accuracy of the E866 experiment is so high that one must 
wonder about the impact of next-to-leading order (NLO) corrections on the 
comparison of Drell-Yan production on different targets. The point is that one 
would like to interpret the Drell-Yan process in terms of the annihilation of 
quarks from the beam hadron on antiquarks in the target hadron. Roughly 
speaking, in a suitable kinematic domain, Drell-Yan production is dominated 
by annihilation of the valence u-quarks of the beam proton on the u-sea in 
the target proton. Because the valence densities are well known, in the pp- 
Drell-Yan production one measures the u sea. Invoking the charge symmetry, 
the pn-Drell-Yan measures the d sea and the ratio of the pn and pp cross 
sections gives an access to the d/u ratio. ^1^0 corrections spoil this simple 
picture. They fall into two broad categoriesH: the QCD radiative corrections 
to the familiar annihilation (ANN) qq ^ J* ^ retain the property of the 
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DrcU-Yan cross section being a (sum over flavours of the) product of quark 
and antiquark densities. The second class of corrections however comes from 
the 'QCD-Compton' (QCDC) subprocesses gq — > (77* and gq — > 97* and has 
the form of a certain convolution of the gluon and quark+ antiquark densities 
in the beam and target and vice versa. The difference in the beam and tar- 
get parton contents of the two classes of corrections and their dependence on 
the form of the parton densities raise several issues: Firstly it is well estab- 
lished that NLO corrections, ususally parametrised in terms of the so-called 
KoY-isiCior.KDY = da dy {LO + N LO) / da dy {LO) , are large: Kdy ~ 1.5 — 2. 
Second, the a;-dependence of the gluon and sea densities is different, and there 
is a troubling possibility that the relative size of the contributions from an- 
nihilation and QCDC NLO corrections changes over the space of the beam 
and target Bjorken variables. Furthermore the NLO corrections depend on 
the shape of parton distributions, and may change the realative importance of 
the uu and dd annihilation contributions compared to the LO formulas. To 
summarize, the issue is that NLO corrections are large, and because in the 
E866 data the statistical accuracy of the measured pn/pp ratio is at the level 
of several percent, a scrutiny of those subtleties of the NLO corrections is called 
uponB. 

The purpose of the present study is an investigation of the beam and 
target dependence of NLO corrections relevant to the interpretation of the 
pn/pp data; we shall be primarily interested in the region of rather large x 
relevant to the E866-data. The contributions from the NLO corrections to the 
differential cross section appear in the form 

^^g- cx ^e^A^j(ii,t2,xi,a;2,Q2) ® Q')/2(t2, Q') + (1 2)] (1) 

for the annihilation corrections of Fig. 1 Ila); and similarly for the QCDC 
corrections (fig 1 IIb,c): 

oc [^fg{hM,x^,x2.Q^)®{^e) [h{ti,Q^) + h{tuQ^)\)g2{h.Q^) 

+ Agf ih ,t2,Xi,X2,Q^)<i^gi ih ,Q^)(J2 4 [/2 (*2 ,Q^) + f2 {t2 , Q')] ) } 

The convolution sign here implies an integration in the variables ti,i = 1,2 
over the domain Xi < U < 1 and all the required coefhcient functions Ajj, A/g 
are available in the literature^. For our purposes it is most convenient to study 
the dependende on flavour etc. in terms of partial /C-factors, that parametrise 
the NLO-corrections in the following manner: 

Affih ,t2,Xl,X2,Q^)® h {h ,Q^)h {t2 , Q') = 
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Afgih,t2,Xi,X2,Q^) ® flitl,Q^)g2it2,Q^) = 

ICfg{xi,nf,X2,ng) ■ fi{xi,Q'^){l - X2)g2{x2,Q^) ■ 

The functions ICj:j,ICfg now bear a dependence on the functional form of the 
parton distributions which is encoded here through the exponents Uf^Ug of 
the large-x behaviour, f{x,Q'^) ~ (1 — x)"f ,g{x,Q^) ~ (1 — x)"". The above 
shown expressions enter the calculation of the experimentally measured cross 
section, and the n / , rig-dependence clearly demonstrates the potential model 
dependence of the extraction of u, d-densities from experimental data. Let us 
emphasize that there is no theoretical a priori reason that the QCDC contri- 
butions should mimick the gg-annihilation processes over a wide (xi, X2, Q^)- 
range. In figure ^ we show a variety of partial /C-factors for two representative 
parametrizations of the parton distributions. A substantial dependence on 
flavour, xi and X2 can clearly be seen. As a matter of fact, for xi,2 > 0.5 — 0.6 
the parton distributions are so poorly known, that the calculations in that 
region have little reliability. In particular the interpretation of Drell-Yan ex- 
periments at larger xi^2 values shall be strongly model dependent. In the region 
of xi^2 relevant to the E866 experiment however the model dependence turns 
out to be weak within the experimental error bars and an extraction of d/u in 
the LO formalism can still be trusted. 
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Figure 1: I: A selection of partial A^-factors calculated with two representative parametriza- 
tions of parton densities; shown as function of X2 for several xi-values: x\ = 0.35 (solid); 
xi = 0.55 (dashed); xi = 0.75 (dot-dashed);i.e. xi grows from bottom to top curves. II: a) 
gluon bremsstrahlung correction to the qq annihilation, b) QCDC contributions. 
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